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A RELIABILITY  GROWTH  MODEL  INVOLVING 
DEPENDENT  COMPONENTS 

by 

N.  Langberg  and  F.  Proschan* 

Abstract 

Earlier  papers  have  shown  how  to  convert  competing  risk  models  involving 
dependent  random  variables  into  models  involving  only  independent  random 
variables,  while  simultaneously  preserving  the  distribution  of  the  minimum  and 
the  probabilities  of  the  various  failure  patterns.  In  the  present  paper,  we 
consider  a sequence  of  such  conversions  occurring  at  successive  points  in 
chronological  time  in  which  the  independent  random  variables  are  becoming 
stochastically  larger.  We  obtain  results  which  essentially  demonstrate  that 
the  limiting  distributions  in  the  sequence  of  dependent  models  "correctly" 
correspond  to  the  limiting  distributions  in  the  sequence  of  independent  models. 

These  results  have  applications  in  reliability  growth  models  and  in  bio- 
medical competing  risk  models  in  which  the  competing  risks  are  increasing 
with  age;  in  these  models  dependency  is  permitted  among  the  random  variables. 
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A RELIABILITY  GROWTH  MODEL  INVOLVING 
DEPENDENT  COMPONENTS 

by 

N.  Langberg  and  F.  Proschan 
Florida  State  University 

1.  Introduction.  In  several  earlier  papers,  Tsiatis  (1975),  Miller  (1977), 
and  Langberg,  Proschan,  and  Quinzi  (1978)  show  how  to  convert  various  conpeting 
risk  models  involving  dependent  random  variables  into  models  involving  only 
independent  random  variables,  while  simultaneously  preserving  the  distribution 
of  the  minimum  and  the  probabilities  corresponding  to  certain  "failure  patterns". 
Explicit  equations  are  presented  which  yield  the  distributions  of  the  independent 
variables.  In  a more  recent  paper,  Langberg,  Proschan,  and  Quinzi  (1977)  develop 
statistical  estimators  of  parameters  of  interest  in  competing  risk  models  in 
which  causes  or  times  of  death  (in  the  biomedical  context)  or  of  failure  (in 
the  reliability  context)  are  not  necessarily  independent. 

In  the  present  paper  we  present  a result  which  should  prove  to  be  basic 
in  converting  reliability  growth  models  involving  dependent  failure  times  into 
equivalent  models  involving  only  independent  random  variables.  Analogous 
applications  exist  in  the  biomedical  field  in  which  survival  functions  may  be 
decreasing,  rather  than  increasing,  with  chronological  time.  The  probabilistic 
theorem  presented  will  be  useful  in  inference  in  reliability  growth  models  involving 
dependent  failure  times. 

^Throughout  the  paper  we  use  "decreasing"  in  place  of  "nonincreasing"  and  "increasine" 
in  place  of  "nondecreasing".  * 


MW— I in'  w”* 


I 


9 


2 


2.  The  Reliability  Grovfth  Model.  The  notation  and  terminology  are  as 
in  Langberg,  Proschan,  and  Quinzi  (1978). 

Consider  a series  system  of  n components  undergoing  improvement  as  time 
passes.  At  fixed  chronological  time  u,  0 s u < uq  (uQ  possibly  infinite) 
component  i has  random  lifelength  T^u)  for  i ■ 1,  .. .,  n,  where  T^u),  ....  Tn(u) 
are  not  necessarily  mutually  independent.  We  say  that  failure  pattern  I occurs 
if  the  simultaneous  failures  of  the  components  in  subset  I of  (1,  . . . , n)  and 
of  no  other  components  causes  (i.e.,  coincides  with)  the  failure  of  the  system. 
Define 


C (100) 


I if  failure  pattern  I occurs  at  time  u 
0 otherwise. 


Let  S(u)  and  T(u)  represent  the  vectors  of  component  life  lengths  of  two  systems 
whose  system  life  lengths  are  S(u)  and  T(u),  respectively.  We  say  that  the 
two  systems  are  equivalent  in  life  lengths  and  patterns  at  time  u 
(S(u)  ==  T(u))  if  P(S(u)  > t,  C(S(u))  - I)  - P(T(u)  > t,  C(T(u))  - I)  for  each 
t * 0 and  each  I c {1,  . ..,  n).  Thus,  two  systems  which  are  equivalent  in  life 
length  and  patterns  are  such  that  (i)  their  life  lengths  have  the  same  distribu- 
tion and  (ii)  the  corresponding  failure  patterns  in  the  two  systems  have  the 
same  probability  of  occurrence. 

The  main  result  of  Langberg,  Proschan,  and  Quinzi  (1978)  may  be  restated 
in  the  context  of  our  model  as  follows: 


2.1.  Theorem.  For  fixed  u,  0 s u < uq,  let  T(u)  - minCT^u),  1 s 1 s n) 

denote  the  life  length  of  an  n-component  series  system,  where  T^ (u)  represents 

the  life  length  of  «nponent  i,  i ■ 1,  n.  Define  F (t)  • P(T(u)  > t, 

u,  i 

C(T(u))  ■ I),  Fu  j(t)  - P(T(u)  s t.  C (T(u) ) - I),  ^(t)  ■ P(T(u)  > t),  and 
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" sup{x:  Fu(x)  > 0>.  Then  the  following  statements  hold: 

(i)  A necessary  and  sufficient  condition  for  the  existence  of  a set  of 
independent  random  variables  (HJ#  I c {l,  . . . , „»  which  satisfy  H(u)  kL  T(u) , 

where  H(u)  ■ min(Hj(u),  I c {1,  . n}),  is  that  the  sets  of  discontinuities 
of  the  F _ be  pairwise  disjoint. 

U 1 1 

(ii)  The  distributions  of  Hj(u),  I « {1,  ...»  n},  are  uniquely  determined 
on  the  interval  [0,  a(Fu))  as  follows: 


Gu,I(t)  " P(Hi<u>  > t) 


•n 

•uf1*^ 


F 

u 


Uu(I,j))/[Fu(au(I,j))  ♦ fu(an(I,j))], 


(2.1) 


° s t -c  a(Fu),  where  F^j  is  the  continuous  part  of  F^j,  (a^I.j)}  is  the  set 

of  discontinuities  of  P^.  in  [0  tl  and  *u  ^(I.j))  is  the  size  of  the  jumps  of  F, 
»t  au(I,j). 

(Note  that  Gu  j may  place  mass  at  infinity.) 

The  survival  probabilities  = l - g j for  Hj(u),  the  time  until  a 
shock  occurs  which  simultaneously  destroys  subset  I,  for  I c {1,  ....  n)  are 
obtained  by  solving  the  identities: 


and 


n 

Ic(l 


»•••! 


(t) 


Fu(t),  t 2 o, 


(2.2) 


{ JII5u,dGu.I  ■ 1 * »• 


(2.5) 
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Thus  at  each  fixed  instant  u of  chronological  time,  we  may  "replace"  the 
original  series  system  of  dependent  components  by  a set  of  mutually  independent 
sources  of  shock;  a given  source  of  shock  fails  a corresponding  subset  of  com- 
ponents. Moreover,  the  replacement  is  so  chosen  as  to  yield  the  same  joint 
distribution  of  failure  pattern  and  system  life  length  as  possessed  by  the 
original  system.  Finally,  the  distribution  of  the  time  of  each  type  of  shock 
occurrence  is  explicitly  specified  in  (2.1). 

Suppose  now  that  the  system  is  experiencing  reliability  growth;  specifically, 
assume  that: 

(i)  G T is  increasing  in  u,  0 £ u < u for  each  I e {l,  ...,  n},  so  that 
Hj(u)  is  stochastically  increasing  in  u,  0 s u < uq  foT  each  I e (l,  . n}. 

That  is,  each  type  of  shock  is  occurring  with  decreasing  frequency  as  chrono- 
logical time  passes.  This,  in  turn,  implies  that  the  system  life  length  T(u) 
is  stochastically  increasing  in  u.  Since  G T(t)  is  monotone  increasing  in 

u < u , it  follows  immediately  that  the  G .(t)  converge  to  inf  sup  G T(x) 

° u’  _ _ x>t  u<u  u* 

for  each  t 2 0,  I c (l,  ....  n};  call  this  limit  Gj(t).  Gj(t)  is  ° a survival 

function,  with  mass  possibly  at  infinity. 

Three  basic  questions  now  arise: 

(1)  Does  the  joint  distribution  F T of  failure  pattern  and  life  length 

u,  1 

converge  to  a joint  distribution  as  u ■*  u#  for  each  I c (1,  ...,  n}? 

(2)  If  such  limiting  joint  distributions  exist,  do  they  coincide  respectively 
with  the  solutions  (call  them  F j)  of  (2.1)  corresponding  to  the  limit  distributions 
Gj,  I c (1,  ....  n)  (as  we  would  hope)? 

(3)  Does  the  system  lifelength  distribution  Fy  converge  to  a lifelength 
distribution  aj  u ♦ u#;  if  so,  does  it  coincide  with  the  solution  F obtained 
from  (2.2)  where  the  subscript  u is  omitted? 

In  the  next  section  we  prove  that  under  mild,  reasonable  conditions,  affirm- 
ative answers  exist  for  all  three  questions. 
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3.  Limit  Results  in  the  Reliability  Growth  Model.  In  this  section  we  give 
convenient  sufficient  conditions  to  yield  affirmative  answers  to  questions  (1), 
(2),  and  (3)  at  the  end  of  Section  2,  and  present  proofs  of  our  results. 

We  are  considering  the  model  of  Section  2.  We  make  the  following  two 
additional  assumptions: 

(ii)  The  sets  of  discontinuities  of  the  G . are  pairwise  disjoint  for  each 

u,  i 

fixed  time  point  u,  0 £ u < uq. 

(iii)  The  sets  of  discontinuities  of  the  Gj  are  pairwise  disjoint. 

The  main  result  may  now  be  stated: 


3.1.  Theorem.  Assume  the  reliability  growth  model  specified  in  Section  2 
and  assume  (i),  (ii),  and  (iii).  Then 


(1)  lim  F j(t)  exists  for  each  t 2 0 and  each  I c (l,  n), 

u-m  ’ 
o 


(2)  lim  F j(t)  * F (t)  for  each  t £ 0 and  each  I c (1,  . ...  n), 

u-m  ' 
o 


(3)  lim  F (t)  exists  and  = F(t)  for  each  t i 0;  moreover,  F (t)  is 

u-*u  u 

o _ 

increasing  upward  to  F(t)  as  u increases  upward  to  uq. 

To  prove  Theorem  3*1,  we  will  make  use  of 


3.2  Lemma . Let  u^  be  an  increasing  sequence  converging  to  uq.  Then  for 
each  t i 0 and  I c (l,  . n),  lim  sup  F T(t)  s F.(t). 

k-*.  V1  1 

Proof  of  Lemma.  By  (2.3),  we  may  write: 


/ n Cf  . dG 

t j.i  vJ 


i . s / n 
V ; j. 


G . dG  _ , 

j.i  J v1 


since  each  G . 

un‘J 


Next  note  that 
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fjl,  G^G“,.i ' - cvi(0>  d*  11  - /„  gj(x)i 

00 

4 l Vw  - \.IW1  d,  » • "nl  'jWi- 

But  by  the  Lebesgue  dominated  convergence  theorem, 

oo 

ii:  l \.'M  - dx  ^ =,(»] 

00 

■ / [Gj(t)  - G(x)]  d [1  - n G.(x)] 

* J*I  J 

00 

- / n G (x)  dG  (x)  * F_(t), 
t J*I  J 1 1 

by  assumption  (iii) . The  desired  conclusion  now  follows.  || 

We  may  now  prove  Theorem  3.1. 

Proo^of  Theorem  3.1  (1),  (2).  Let  (u^,  0 s u^  < uq,  be  an  arbitrary 
sequence  converging  to  uq  and  let  t i 0.  By  (2.2), 


I 


F..  TCt)  * n 


, l TV«-/  - II  G _(t)  * 

*<1 n)  V1  Ic{l, . . . ,n)  uk’1^ 


But 


lim  sup  n G _(t)  » lim  sup  V F (t) 

k~  n}  V1  ^ Ic(i,L,n>  VI™ 


sup  T ,,  l ,FU  !<*)  s l lim  sup  F T(t)  s T p 
ic{i,...,n>  uk*J  mi n>k~>  uk,i  Ic{i,:..,n} 


by  Lenma  3.2.  Now  by  assumption  (iii)  and  the  appropriate  version  of  (2.2) 
have: 

l Fr(t)  - n frT(t). 

Ic{l,...,n)  Ic{l,...,n)  1 


we 


j(t). 
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Recalling  that  Gj(t)  » lim  j(t:)  from  Section  2 and  summarizing  the  inequalities 
above,  we  may  state: 


Ic(l,...,n} 


GT(t)  S 


l 


lim 

,n}k-*« 


sup 


V(,) 


l 

I<=(l,...,n} 


Fj(t) 


n 

lc(l,. 


,n} 


GjCt) 


The  desired  conclusions  (1)  and  (2)  now  follow  immediately. 


Proof  of  Theorem  3,1  (3). 


exists.  Next  note  that 


By  (2.2),  lim  F (t)  exists 
u-m 

o 


since  each  lim 
u-m 

o 


F(t)  = T I ,FI(t)  = n (t)  * lim  It 

Ic(l,...,n)  u-m  Ic{l 

A * ' 


. ,n) 


lim 

u-m 


l F T(t)  = lim  F (t). 
Ic( 1 , . . . ,n)  u* 1 u^i  u 


The  monotonicity  of  Fu(t)  in  u is  a consequence  of  the  monotonicity  of  each 


*»,r 
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4.  Extensions,  Modifications,  and  Generalizations. 

4.1.  Remark . In  the  reliability  growth  model  discussed  above,  we  assumed 
that  reliability  growth  was  occurring  during  the  interval  of  chronological 
time  [0,  tQ) . Obviously,  the  conclusions  of  Theorem  3.1  may  be  obtained  by 
assuming  that  reliability  growth  occurs  at  the  time  points  of  any  subset  U 

of  [0,  ■»)  (or,  for  that  matter,  over  the  entire  half-line).  As  examples, 
consider: 

(a)  U = (u,  2u,  3u,  where  u > 0, 

(b)  U = (0  s Uj  < < u3  < ...} 

(c)  U <=  { [0,  uq]  < u1  < u2  < . . .} , 

etc.  Moreover,  the  set  U need  not  be  deterministic,  but  may  be  random. 

4.2.  Remark.  In  the  model  of  Section  2 and  the  corresponding  result. 

Theorem  3.1,  we  assumed  that  shock  intervals  from  each  source  were  stochastically 
increasing  with  chronological  time.  It  is  clear  that  by  reversing  inequalities 
and  the  direction  of  monotonicity,  we  can  obtain  a dual  to  Theorem  3.1  in  which 
the  intervals  between  shocks  from  each  source  are  stochastically  decreasing 

with  chronological  time.  System  reliability  (or  in  the  biomedical  context, 
organism  survival  probability)  would  decrease  with  chronological  time. 

4.3.  Remark.  The  model  of  Section  2 and  the  limit  results  of  Section  3 
were  formulated  in  terms  of  reliability  growth.  An  equally  important  and  useful 
application  exists  in  the  context  of  competing  risks  of  death  of  a biological 
organism  from  a variety  of  diseases,  accidents,  or  other  causes.  In  this 
context.  Remark  4.2  is  especially  relevant.  We  consider  a situation  in  which 

an  organism  is  subject  to  death  as  a result  of  any  of  a number  of  diseases 
or  combinations  of  diseases.  As  time  passes,  the  organism  becomes  more 
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vulnerable  to  each  of  the  diseases.  Note  that  the  model  now  appropriate  differs 
from  the  model  of  Section  2 only  in  the  direction  of  the  monotonicity;  i.e., 
the  different  types  of  shocks  are  becoming  stochastically  more  frequent,  rather 
than,  as  in  the  first  model,  stochastically  less  frequent. 

4.4.  Remark.  In  the  reliability  growth  model,  we  are  able  to  deduce  mono- 

tonocity  of  Fu,  system  reliability  as  a function  of  chronological  time  u. 

However  it  is  not^  necessarily  true  that  each  F . is  monotone  in  u.  This  is  a 

u , x 

consequence  of  the  fact  that  as  some  modes  of  failure  become  less  likely,  others 
may  become  more  likely. 

4.5.  Remark . Assumption  (iii)  of  Theorem  3.1  can  not  be  obtained  as  a 
consequence  of  assumptions  (i)  and  (ii) . A simple  counterexample  is  available 
to  verify  this  assertion. 
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W.  ABSTRACT 

Earlier  papers  have  shown  how  to  convert  competing  risk  models  involving  dependent 
random  variables  into  models  involving  only  independent  random  variables,  while  simul- 
taneously  preserving  the  distribution  of  the  minimum  and  the  probabilities  of  the  various 
failure  patterns.  In  the  present  paper,  we  consider  a sequence  of  such  conversions  occurring 
at  successive  points  in  chronological  time  in  which  the  independent  random  variables  are 
becoming  stochastically  larger.  We  obtain  results  which  essentially  demonstrate  that  the 
limiting  distributions  in  the  sequence  of  dependent  models  "correctly"  correspons  to  the 
limiting  distributions  in  the  sequence  of  independent  models. 

These  results  have  applications  in  reliability  growth  models  and  in  biomedical  competing 
risk  models  in  which  the  competing  risks  are  increasing  with  age;  in  these  models  dependency 
is  permitted  among  the  random  variables. 


